Direct demonstration of circulating currents in a controllable 7T-SQUID generated by 

a to 7r transition of the weak links. 
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A controllable 7T-SQUID is a DC SQUID with two controllable 7r-junctions as weak links. A 
controllable -n -junction consists of a superconducting - normal metal - superconducting Josephson 
junction with two additional contacts to the normal region of the junction. By applying a voltage V c 
over these contacts it is possible to control the sate of the junction, i.e. a conventional (0) state or 
a 7r-state, depending on the magnitude of V c . We demonstrate experimentally that, by putting one 
junction into a 7r-state, a screening current is generated around the SQUID loop at integer external 
flux. To be able to do this, we have fabricated controllable 7r-junctions, based on Cu-Nb or Ag-Nb, 
in a new geometry. We show that at 1.4 K only the Nb-Ag device shows the transition to a 7r-state 
as a function of V c consistent with theoretical predictions. In a controllable n SQUID based on 
Nb-Ag we observe, a part from a screening current at integer external flux, a phase shift of n of the 
Vsquid — B oscillations under suitable current bias, depending on the magnitude of V c - 
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I. INTRODUCTION 

In recent years considerable attention has been paid 
to controllable Josephson junctions. In general, such a 
Josephson junction consists of a Superconductor-normal 
metal-superconductor (SNS) junction with additional 
electrodes connected to the normal region of thew junc- 
tion. Inside the normal region coherent quantum states 
are formed due to the presence of the two supercon- 
ducting electrodes. The properties of the junction, and 
in particular the critical current I c are determined by 
the energy spectrum and the occupation of these states. 
Hence it is possible to change I c by changing the energy 
distrih^k»|-pf the quasiparticles in the junction normal 
rcgionBoaula. This is done by sending a control current 
through the additional electrodes connected to the nor- 
mal region or by applying a control voltage V c over them. 
Several experiments have been performed the past years 
on these devices, which differ mainly in thepWayJ^which 
the energy of the quasiparticles is changedOLrljoEj. The 
most interesting situation is that of the so called con- 
trollable 7r-junctionsBEj. In these systems it is not only 
possible to change the magnitude of the critical current of 
the junction but also to reverse its direction with respect 
to the phase difference ip between the two superconduct- 
ing electrodes. This sign reversal is equivalent to the 
introduction of an extra phase factor tt in the Josephson 
supercurrent (J sc )-phase relation 

I sc = I c sin(tp) => I sc = -I c sin{ip) = I c sin(<p + it) . 

(1) 

This means that for zero supercurrent, the ground 
state is that of a macroscopic phase difference equal 



to 7r. There is an analogy here with 7r-junctions 
associated with bi-crystalsEJ or 's-d' contacts in ce- 
ramic superconductors^ and 7r-junctions using a dilute 
ferromasp.et| f as the "normal" region of a S-f'-S 
junctionI3E-X However, the ceramic systems remain a 
7r-junction (or normal junction) once fabricated and the 
S-f'-S junctions switch between a normal and a 7r-state 
only as a function of temperature. This in contrast 
with the controllable 7r-junction where the state of the 
junction as well as the magnitude of the critical current 
is fully controllable by means of the control voltage. 

Using such controllable 7r-junctions as weak links 
of a DC SQUID (direct current superconducting quan- 
tum interference device) , it is possible to tune the critical 
currents of the individual weak links. More interestingly, 
it is possible to add an extra phase of n or 2tt to the 
SQUID loop without the application of an external 
magnetic field by switching one or two of the junctions 
in a 7r state. An extra phase factor of ir is identical to 
the external application of half a flux quantum, 4>o/2, 
with (j>o — 2.07 • 1CU 7 Wb. This is a direct consequence 
of the condition of a single- valued wave function around 
the SQUID loop, 



2nc 



— ipi — ip2 = 2im , 



(2) 



where <fi is the total flux in the ring, n is an integer 
and ifi and if 2 are the phases over the two Josephson 
junctions, as shown in Fig.|l]a. Switching one junction in 
the 7r-state at an external flux <fi = n(f>o adds tt to the left 
hand side of Eq||. This leads to the same solution of the 
equation if the junction would be in the "conventional" 



state with an external flux 4> — [n ± -k)4>o- Hence it 
is possible to generate a screening current in a SQUID 
with a sufficiently large inductance in the absence of a 
magnetic field. 

The main purpose of this paper is to demonstrate 
unequivocally that it is possible to add a macroscopic 
quantum phase of it or 2ir to a DC SQUID with 2 
controllable 7r-junctions as weak links. We demonstrate 
that this results in a screening current that can be 
switched on or off depending on V c at an external flux of 
an integer flux quantum. 

The outline is the following: We begin, in section 
II, with a very short description of the basics of the DC 
SQUID and a description of the theoretical concepts 
regarding controllable 7r-junctions. In section III we dis- 
cuss the sample design and fabrication and in section IV 
we discuss experiments on single controllable 7r-junctions 
made of Nb-Ag or Nb-Cu. The geometry of these devices 
differs from the conventional cross shape. We show that 
the Nb-Ag device shows a to ir transition at 1.4 K, in 
contrast with the Nb-Cu device. The full dependence 
of the critical currents of both devices on V c is in good 
quantitative agreement with the theoretical predictions 
discussed in section II. In section V we proceed with 
measurements on a controllable SQUID, made of Ag-Nb. 
We show that the magnitude and sign of the voltage 
oscillations of the SQUID as a function of the magnetic 
field under suitable current bias can be tuned. In the 
last section we measure directly the flux expelled by the 
screening current in a controllable SQUID by placing a 
smaller SQUID directly behind it. We show that the 
screening current around integer external flux can be 
switched on or off by switching one of the junctions of 
the controllable SQUID in a tt- or a normal state. 



II. THEORETICAL CONCEPTS 
A. the DC SQUID 

Before we discuss the controllable 7r-junction.it.. ma 
helpful to discuss the basics of the DC SQUIDBMI 1 



be 
A 

DC SQUID is a ring of superconducting material with self 
inductance L and two Josephson junctions with a critical 
currents l c \ and l c i . In the absence of any external flux, 
(f) e xt, the SQUID is in equilibrium and behaves as a sin- 
gle Josephson junction with a total critical current I c = 
Ici+Ic2- If 4>ext 7^ 0, then the total flux in the SQUID 
loop changes, depending on its self inductance, according 
to 



LJ 



(3) 



with J the screening current running around the SQUID 
loop. The flux in the loop, <p, also depends on the phases 
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FIG. 1: Basic operation of the DC SQUID. (A) shows a DC 
SUID with a screening current J due to an external flux (j) ex t- 
The supercurrent of the SQUID, as shown in (B), oscillates 
as a function of cj> e xt- By applying a current bias Ib the 
oscillating critical current is transferred into an oscillating 
voltage Vsquid between contacts 1 and 2 (A), as shown in 
(D). The screening current J has, in the limit of large self 
inductance L, a saw tooth shape as a function of <f>o, as shown 
in (C). 



of the two Josephson junctions 

2tT0 
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(4) 



which is directly obtained from Eq.^[ In the limit of small 
self inductance, so that the LJ term in EqJ^ is negligible 
and if I c \ = I c \ = I c i,2, it is easy to show that 
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where we use Eqs.[j]]| an d ||. In this case I c is periodic 
in (jj and varies between 2/ cl 2 at c/) ext — n<f>o and at 
4>ext = (n + 1/2)00- At nonzero inductance the period- 
icity remains the same but the amplitude of the I c mod- 
ulation is reduced, for I c —0 is not reached at (j> = it/2 
due to the presence of the circulating current J. In the 
limit of large self inductance, the amplitude of the I c 
modulation decreases to <fio/L, whereas J becomes a saw 
tooth function varying linearly from —<po/2L to +0o/2L, 
as shown in Fig.|l|(B) and (C). The periodic behavior of 
I c as function of <p ex t can be transformed into a periodic 
voltage by sending a bias current Ib a little larger than 
the maximum value of I c through the SQUID. Due to 
the non-linear I-V curve of the SQUID the difference in 
Ic{<t>ext) translates into a voltage VsQuiD^ext), thus ere- 
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FIG. 2: Schematic drawing of a controllable 7r-junction. The 
center of the control channel, which connects two large reser- 
voirs, is coupled to the normal region of a co-planar SNS 
junction. By applying a voltage V c over the reservoirs the 
energy distribution of the quasiparticles in the channel and 
hence in the normal region of the SNS junction is changed. 
As a result the occupation of the suercurrent carrying states 
in N is changed, resulting in a change in the critical value of 
the supercurrent Isns- 



ating a very sensitive flux to voltage transformer. The 
resulting Vsquid — B oscillations are shown in Fig.[T](D). 



B. Controllable 7r-junctions 

A controllable 7r-junctions is a device that combines 
two different structures in one: A thin film SNS junction 
and a short mesoscopic wire (which we call control chan- 
nel) between two large electron reservoirs in thermal 
equilibrium. A schematic drawing is given in Fig.|: The 
center of the control channel, which connects two large 
reservoirs, is coupled to the normal region of a co-planar 
thin film SNS junction, in which the normal region 
partly overlaps the superconducting electrodes. The 
length and width of the normal region are denoted by L 
and W in the figure. The size of the junction L < 1/xm 
and the control channel length, I ~ 5/im, is larger than 
the elastic mean free path of the electrons, which is 
approximately given by the film thickness (~ 50nm). 
Hence the electron motion is diffusive. 

The principle of operation is the following: The 
presence of the two superconducting electrodes induces 
a continuous spectrum of supercurrent carrying states 
Im(J(E)) in the normal region of the junction, which 
is responsible for the supercurrent transport (Isns in 
Fig.||). The energy scale E is the energy with respect 
to the Fermi level. As stated before, the critical value 
of Isns, Ic, is determined by the energy spectrum and 
the occupation of these supercurrent carrying states. 
The occupation of the states is determined by the 
energy distribution of the quasiparticles in the normal 



region of the junction, f(E), which is identical to the 
energy distribution in the center of the control channel. 
This energy distribution can be modified by means 
of applying a voltage V c over the control channel, so 
f(E) = f(E,V c ). The^ela^n between l c , Im(J(E)) 
and f(E, V c ) is given bjffit 



/c = 7T /" dE[l ~ 2f{E > V c)] Im (J(E, V = tt/2)) , 

Jtn J — oo 



(6) 



where (p the phase difference between the superconduct- 
ing electrodesEEl. ft is important to realize that the crit- 
ical current is modified by changing the energy of the 
quasiparticles in the normal region and not by the flow 
of the control current associated with V c . In the present 
geometry this is quite clear, because the current paths 
of Isns and the control current are separated. This is 
in contrast with the conventional cross geometry, where 
Isns crosses the control current. However, even in this 
case it is the energy of the quasiparticles and not the 
control current that determines the critical current. 

To calculate the dependence of I c as a function of V c 
we need to know both the shape of the supercurrent 
carrying density of states, Im(J(E)) and the shape of 
the quasiparticle energy distribution in the normal re- 
gion of the junction, f(E, V c ). The exact energy depen- 
dence of Im(J(E)) is calculated using the quasi-classical 
Green's function theorju-oBa. ft depends on the ra- 
tio of the Thouless energy with the energy gap A of 
the superconductors. The Thouless energy is given by 
E t h — HD/L 2 , withpL the length of the normal region 
of the SNS junctiontj and D the diffusion constant of 
the normal metal. The general shape of Im(J(E)), as 
shown by the solid black line in both panels of Fig. || 
is a strongly damped oscillation with a hard gap at low 
energies E < E t h- The positive and negative parts of 
the supercurrent carrying density of states represent en- 
ergy dependent contributions to the supercurrent in the 
positive and negative direction. 

The other relevant quantity of the system is the quasi- 
particle energy distribution, f(E,V c ). The exact shape 
of f(E,V c ) in the center of the control channel (and 
hence in the junction normal region) depends, as shown 
by Pothicr et a/B3, on the voltage V c applied over the 
control channel and on the amount of electron-electron 
and electron-phonon scattering. At low temperatures, 
T < IK, electron-electron scattering is the dominant re- 
laxation mechanism and the electron energy distribution 
depends only on V c and the electron-electron interaction. 
In this case, we can distinguish two limiting possibilities: 

1: no electron- electron interaction If the control 
channel is sufficiently short, so that the electron-electron 
interaction time To exceeds the diffusion time td through 
the channel, the electrons conserve their energy. The 
resulting electron energy distribution in the center of 
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FIG. 3: Supercurrent carrying density of states Im(J(E)) for 
A/E th =70 at ip = tt/2 (solid line) together with 1 - 2f(E) 
(dotted line) in the case of a thermal distribution (a), and 
a step distribution (b). The shaded area represents the con- 
tributions to the integrand of Eq.|(| E represents the energy 
with respect to the Fermi level an Eth is the Thouless energy. 



the wire is given b; 



f(E,V c ) = 1/2 



fd[E 



2 



fd\E 



2 



(7) 



here fd(E) is the Fermi-Dirac distribution at a tempera- 
ture T 



fd(E) = 



exp 



E 



(8) 



with Ub Boltzmann's constant. Eq.^ is the renormalized 
addition of the thermal distribution functions of the 
reservoirs with an energy separation of eV c . This 
represents, if khT < eV c , a double step structure with 
an occupation fraction of 1/2 in an energy range of eV c 
around the Fermi energy. Such a step distribution will, 
at a large enough value of V c , block all the positive 
contributions of Im(J(E)), as shown in Fig. ||b, because 
an occupation fraction of 1/2 causes the integrand of 
Eq.|| to vanish. This results in a reversal in the direction 
of the supercurrent and hence the transition to a 7r-state 

2: very strong electron- electron interaction If the control 
channel is very long, so that To <C 173, the electron 
system regains a local thermal equilibrium. The energy 
distribution given by the Fermi-Dirac distribution, with 
an effective temperature T e ff > Tsathi 



due to the compensation of the positive and negative 
contributions in Im(J(E)). 

The shape of the distribution function in a real 
experiment is in general in between these two limits, 
with possible modifications due to elect ron-phonon 
interactions. We calculate the exact shape of the energy 
distribution function using a model presented by Pothier 
et alEjlEl. In this model the shape of the electron energy 
distribution in the wire depends on four parameters: 
The voltage V c applied over the wire, the effective 
strength of the electron - electron interaction, the 
effective strength of the electron phonon interaction and 
possible heating of the reservoirs. The electron-electron 
interaction is described by a two-particle interaction 
with an interaction Kernel K a e~ a ■ to- Here e is the 
energy transferred between two interacting particles and 
td the diffusion time through the wire. The magnitude 
K a and energy dependence a of the electron - electron 
interaction can be calculated by means of a direct 
calculation of the screened Coulomb interaction .iu a 
homogenous 1-dimensional diffusive conductorE-SEl 
Without going into detail it is important to realize that 
a direct measurement of the electron energy distribution 
function in diffusive ID gold, copper.. and, silver wires us- 
ing tunnel junction spectroscopyESEjtJ'LJ yields material 
dependent results for K a e~ a which are not_pxmsistent 
with the theoretical predictions using refs.Ej'LJ alone. 
Because the control channels used in oju^pqjieriment are 
similar to the wires measured in refs-EHttlrtJO, we have 
used the experimental values of a in our calculations: 
a = 2 for copper and a = 1.2 for Ag. The strength of 
the interaction, K a is used as a fit parameter. 

In the model, the effect of electron-phonon scattering 
is implemented in a similar way. The magnitude of the 
interaction kernel is, in this case, deduced from mea- 
surements of the phase relaxation timoEa, and given by 
Kphonon — Sns~ 1 meV~ 3 for Cu and Ag. 

As a last input parameter, possible heating of the reser- 
voirs should be taken into account. The power associated 
with the application of the control voltage, P = V c 2 /i? c 
(R c is the control channel resistance) is injected into the 
reservoirs. At low bath temperatures (T < IK), low val- 
ues of R c and large values of V c > lmV, the electron 
gas in the reservoirs can attain a much . h i gh er effective 
temperature than the bath temperatureEZE3. However, 
all experiments presented in the remainder of the text 
are performed at relatively high temperatures_cjf 1.4 K, 
resulting in negligible heating of the reservoirso. 



T, 



eff 



T, 



bath 



+ {aV Ci 



(9) 



Here Tb a th is the bath temperature and a can be 
calculated using the Wiedemann-Franz law, resulting in 
a=3.2 K/mV. This situation is called the 'hot electron 
limit'. The result, a Fermi-Dirac distribution with 
a width increasing with V c , causes a rapid decay of 
the supercurrent, as shown in the left panel of Fig.0a 



III. SAMPLE DESIGN AND FABRICATION 

So farp-pall experimental evidence of a controllable 7r- 
junctioroEj has been achieved using devices with a cross 
geometry. The control channel and the normal region 
of the junction are in fact a single cross shaped metal- 
lic wire, with superconducting contacts at the first two 
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FIG. 4: Scanning electron micrograph of a controllable n- 
junction using side contacts. In this case the normal metal is 
silver, whereas also similar samples using copper are made. 



opposite ends and large reservoirs at the others. This 
geometry is not suitable for use in a SQUID geometry, 
because one of the reservoirs of each junction would have 
to be within the SQUID loop, which is not possible due 
to the large reservoir sizes needed. Therefore we use a 
device geometry with side contacts (see Fig.|j), such as 
used by Morpurgo et al.a in the thermal limit. The dis- 
advantage of this geometry is that it imposes a larger 
minimum length on the control channel than the conven- 
tional cross geometry, where this length can be reduced 
to about 1 /im. As a consequence, a material with a 
slow electron-electron relaxation is needed to be able to 
maintain a non-thermal energy distribution in the con- 
trol channel. . For , this , r eason gold is unsuitable as the 
normal metalBOBBEl and copper or preferably silver 
should be used. 

We now describe the sample fabrication, referring to 
Fig. ^. The samples have been realized on a thermally 
oxidized Si wafer that is covered with a 150 nm layer of 
sputter deposited AI2O3 to improve the adhesion of Ag 
and Cu. In the first step 50 nm of Nb is deposited using 
standard e-beam lithography, DC sputtering and subse- 
quent lift-off. The critical temperature of the sputtered 
film is 8.1 K. Subsequently the normal region, control 
channel and the reservoirs (Ag or Cu) are deposited using 
shadow evaporation without breaking the vacuum. This 
is needed because the adhesion of these metals is poor, 
implying that it is not possible to bake a deposited film to 
be able to do another lithography step. We use a double 
layer of PMMA-MA and PMMA with e-beam lithogra- 
phy and wet etching to create a PMMA suspended mask. 



TABLE I: Parameters of the single SNS junctions. Rn is the 
normal state resistance of the SNS junction, L is the length 
of the normal region and W the width (see Fig.^) . E t h is the 
Thouless energy obtained using D and L. 



sample 


L x W [nm] 


R„ M 


D [m 2 /s] 


E th [»eV] 


Ag 


1000x400 


0.65 


0.023 


14 


Cu 


800 x 600 


0.8 


0.021 


21 



The deposition of the metals is done in an UHV deposi- 
tion system with a background pressure of 5 TO -10 mBar, 
the pressure in the system during the evaporation steps is 
< 5 • 10~ 8 mBar. Prior to deposition we use Argon etch- 
ing {Pat — 1 ' 10~ 4 mBar, 500 V) for 3.5 minutes to clean 
the Nb surface. Subsequently we deposit, in the case of 
an Ag sample, 10 nm of Ti adhesion layer under a large 
angle (47°, resulting in a shift to the left in the figure), 
with the result that the Ti layer is only deposited on the 
substrate at the position of the reservoirs and not at the 
position of the thin openings defining the control channel 
and the normal region of the junction. This is important, 
for the effect of a thin Ti layer upon the electron-electron 
interaction time is not known. Then we deposit 50 nm 
Ag perpendicular to the substrate, thus creating the con- 
trol channel and the normal region of the junction. As 
a last step 700 nm of Ag is deposited again at 47° to 
form the reservoirs with an effective thickness of 475 nm. 
In Fig.[| the left shift due to the large evaporation an- 
gle of the reservoirs with respect to the control channel 
is clearly visible. A (parasitic) projection of the control 
channel and the normal region of the SNS junction is not 
present, because the aspect ratio of the top resist is so 
large for these structures that the metal is deposited into 
the side wall of the resist. This results in a removal of 
these projections during lift-off. The fabrication of the 
Cu samples is identical except for the fact that Cu is used 
for the normal conducting regions of the device and that 
no Ti adhesion layer is used underneath the reservoirs. 
We found that this is not necessary for sufficient adhe- 
sion. The length of the control channel cannot be made 
shorter, because the resist will become unstable if the 
large openings defining the reservoirs come any closertil. 
To measure the quality of the niobium - normal metal in- 
terface we have made, in every fabrication run, a cross of 
a 200 nm wide Nb and normal metal wire. The resistance 
of the 200x200 nm interface has been determined to be 
ps 0.1 Q each time, which is smaller than the square re- 
sistance of the normal metal (0.4 f2), indicating that the 
interface is clean. 



IV. EXPERIMENTS ON SINGLE JUNCTIONS 

In the experiment, performed at 1.4 K in a pumped He 
bath, we measure the critical current of the SNS junction 
(from contacts 1 to 2) as a function of the voltage V c 
applied over the control channel (see Fig. 0). The mea- 



6 



TABLE II: Result of the fits on the single-ijitwnctions. For a 
we use the experimental values from Ref.E9Ej: q=2 for Cu 
and 1.2 for silver. 



sample 


K a [ns- 1 fieV a - 2 ] 


E th [neV\ 


td [ns] 


Ag 


0.65 


13 


2.3 


Cu 


1.7 


21.5 


2.5 



surement of the critical current consists of sweeping the 
bias current while measuring the differential resistance 
dV/dl using a small AC modulation on the bias cur- 
rent. We define the experimental critical current where 
the differential resistance reaches half its normal state 
valudlj. The control voltage is applied using another 
current source, with floating ground and V c is measured 
by means of voltage probes close to the control chan- 
nel ends (not visible in the figure). The measurement of 
V c shows a slight modulation if the current bias through 
the SNS junction is smaller than the G-ritieal current, as- 
sociated with Andreev interferometryL3E2l. This modu- 
lation is used to determine the state of the junctiona. 
The experiment is performed on several devices made of 
Cu-Nb and Ag-Nb. All samples show identical behavior 
and we present the data of two samples, one made of 
Ag-Nb and one made of Cu-Nb. The exact parameters 
of these two samples are given in Table |. The experi- 
mental critical current as a function of V c is shown by 
the solid circles in Fig. ^. From the figure is is obvious 
that the Ag-Nb device shows a transition to a 7r-junction 
above V c = V c ,criticai — 430^y and that the device made 
of Cu-Nb does not show this transition. The magni- 
tude of the equilibrium critical current at 1.4 K is, in 
the case of Ag, in good agreement theoretical predictions 
presented by Dubos et alrB. The measured fraction of 
I c R n /E t h — 0.44, whereas 0.5 is predicted. For the cop- 
per sample the agreement is less good. In this case we 
measure 0.5 whereas 2 is predicted. The next step in our 
analysis is a quantitative comparison of the experimental 
data with the theory. As input parameters to calculate 
the distribution function,jye-juse, for a and K p h onons the 
values given in section IlElEil, K p honons=& ns~ 1 meV~~ 3 , 
a=1.2 for Ag and 2 for Cu. As stated before, reservoir 
heating is negligible due to the relatively high bath tem- 
perature and large control channel resistances Xs^8 O for 
both samples), despite the large values of V c fB. Using 
E t h and K a i p h a as fit parameters, we obtain the best fit, 
as shown by the black lines in both panels of Fig.[|, using 
the parameters shown in Table ||. In the Table, we have 
calculated T£> = using for 1 the total diffusion length, 
given by the length of the v-shaped wire plus two times 
the T-shaped extension towards the junction, yielding 
1=5+2- l=7//m (see Fig.|]. The fits are very good over 
the entire energy range. Moreover, the results on K a are 
in good agree ment . with the values obtained in the ex- 
periments fromEIEil, in contrast to pr-pmilar analysis on 
a conventional 7r-junction using golcEJ. The values of the 
Thouless energy inferred from the fits are both in good 



agreement with the geometrical values. The solid lines in 
the inserts of the figure show the electron distributions, 
inferred from the fits, at V c =0.8 mV. For both samples 
the distributions are relatively rounded. In the case of 
the Ag device, this is mainly caused by the bath temper- 
ature of 1.4 K. In the case of the Cu sample the stronger 
electron-electron interaction is responsible for the extra 
rounding. However, the distribution is, even in this case, 
not quite a thermal one but still has a weak double step 
signature. This can be seen by the difference between 
the solid and the dotted line, which represents a ther- 
mal distribution with T e ff — 2.9K. This temperature 
corresponds to the effective temperature in the hot elec- 
tron regime, using Eqj9|with U c =0.8 mV and T( )ot / l =1.4K. 
Model calculations indicate that, in the case of the Cu- 
Nb device, a slight decrease in control channel length or 
measuring at a lower bath temperature would result in a 
distribution that is sufficiently less rounded. This would 
resuLtjn a transition to a 7r-state at large enough values 

of vM 

We conclude that, in the present geometry, only a Nb- 
Ag junction shows the transition to a 7r-junction at 1.4 
K. The absolute value of the critical current and the re- 
sponse of the critical current to V c are all consistent with 
the experimental and theoretical data available. A Nb- 
Cu junction does not show a transition to a 7r-junction 
at this temperature, but would probably do so at lower 
temperatures. To be able to measure at 1.4 K and higher 
temperatures we have fabricated controllable 7T-SQUIDS, 
DC SQUIDS with controllable n- junctions as the weak 
links, from Nb-Ag. 



V. CONTROLLABLE tt-SQUID 

A practical realization of a controllable 7T-SQUID made 
of Nb-Ag is shown in Fig. |^. The fabrication procedure 
is identical to the fabrication of the single junctions de- 
scribed in Section III. The length L of the normal regions 
of both junctions is 1100 nm and the width W of the nor- 
mal regions is 520 nm for the top junction and 220 nm for 
the bottom junction. The resistance of the control chan- 
nels is 8 O, yielding a diffusion constant D=0.023 m 2 / s. 
From L and D, E t h is estimated to be 12.6 neV for both 
junctions. The area of the SQUID loop is 12/im 2 . The 
Thouless energy of both junctions of the SQUID are al- 
most identical to the Thouless energy of the single junc- 
tion made of silver discussed in the previous paragraph. 
We therefore conclude that the response of the critical 
current of each junction to V c is comparable to the one 
presented in the top panel of Fig. |[ 



A. Experiments on a controllable SQUID 

In the first experiment we measure the Vsqjjid ~ B 
oscillations of the controllable SQUID at 1.4 K. We 
bias the SQUID with a low frequency AC bias current 
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FIG. 5: Critical current of a controllable 7r-junction as shown 
in Fig.§, made from Cu-Nb (bottom) and Ag-Nb (top). The 
dots represent the data and the solid line the calculated form 
using Eq. ^[ The solid lines in the inset show the electron dis- 
tribution functions at 14=0.8 mV used in the fit. For compar- 
ison, the Fermi-Dirac distribution in the hot electron regime, 
with T e ff — 2.9K, is given by the dotted line. 



{iBias, 80Hz, 4fiA) between contacts 1-2. The ampli- 
tude of the bias current is a little larger than the max- 
imum critical current of the SQUID. The voltage over 
the SQUID is measured, Vsquid, from contacts 1 to 
2, as a function of the applied magnetic field B using 
a lock-in amplifier. This lock-in technique strongly re- 
duces the noise compared to the conventional measure- 
ment method using a DC bias. Simultaneously we send 
a DC current through the bottom control channel and 
measure the resulting control voltage V C 2, using current 
sources with floating grounds. The result, using the de- 
vice shown in Fig is shown in Fig |^. The solid lines 
represents the Vsquid-B oscillations for increasing val- 
ues of V C 2 (Vd = 0). At first the amplitude of the oscil- 
lations decreases with increasing V C 2 and reaches zero at 
U C 2. C ritica/=0.48 mV, indicating that the critical current 
of the bottom junction is equal to 0. At higher values 
of V c .2 the Vsqu i d-B oscillations re-appear, with a shift 
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FIG. 6: Scanning electron micrograph picture of a controllable 
7T-SQUID. 
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FIG. 7: The voltage over the controllable SQUID as a func- 
tion of the external magnetic field B using a current bias 
Ib — 4/j.A (curves offset for clarity). Solid lines: Vd—O, 
U c2 =0.30,0.38,0.46,0.48,0. 54,0.70,2. 7 mV (bottom to top). 



7r in phase with respect to the oscillations at lower val- 
ues of V C 2- The bottom junction and hence the SQUID, 
are now in the 7r state. At zero field we now measure a 
voltage maximum in stead of a minimum. 

Using an additional current source it is possible to ap- 
ply a control voltage also to the top junction. Now it 
is possible to add a phase of 0, n or 2-7T to the SQUID 
loop, as shown in Fig.||. The solid black line represent 
the 0-0 state, where both V c i and V C 2 are smaller than 
Vc,criticai- By putting junction 1 (dotted black line in 
FigM) or junction 2 (solid grey line in the figure) in a 
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FIG. 8: The voltage over the controllable SQUID as a func- 
tion of the external magnetic field B (curves offset for clarity) 
using both controllable junctions. The solid black line rep- 
resents a 0-0 state, the solid grey line and the dotted black 
line represent the situation with only junction 2 or 1 in the n- 
state, respectively. The dotted grey line gives the situations 
for both junctions in a n state. Note that it has the same 
phase as the solid black line. 

7r-state by applying a sufficiently large value of V c we ob- 
serve in both cases an identical shift of n in phase. By 
putting both junctions in a n state, as shown by the dot- 
ted grey line in FigJ|, we regain the original phase of the 
oscillations due to the 2tt periodicity of the macroscopic 
wave function around the SQUID loop. Similar measure- 
ments at 100 mK in a dilution refrigerator and at 4.2 K 
have shown similar resultsa, with however different am- 
plitudes of the V-B oscillations due to the temperature 
dependence of the critical current of the Josephson junc- 
tions. These results indicate that behavior of the SQUID 
in the 7r-state at 4>ext = nfio is identical to that of a con- 
ventional SQUID at (p ext — (n + 1/2) fa- This should 
also be valid for the screening currents in a SQUID with 
a sufficiently large self inductance. 



VI. DIRECT DETECTION OF THE 
CIRCULATING CURRENT IN A 
CONTROLLABLE tt-SQUID 

In this Section we describe an experiment to detect 
the screening current in a controllable SQUID directly 
by placing a conventional SQUID with small self induc- 
tance inside the loop of the controllable SQUID. The de- 
vice is shown in Fig.^j. A controllable SQUID with a 
large area (70^tm 2 ) is placed upon a conventional SQUID, 
which we call the measurement SQUID, with an area of 
15 /iTO 2 . This SQUID has conventional SNS junctions 
made of gold as weak links. The two devices are elec- 
trically isolated from each other by means of 170 nm of 
AI2O3. It is important to realize that the design of the 
two SQUIDS is critical: The product of the screening 



current J and the self inductance L of the measurement 
SQUID should be negligible, to prevent a disturbance of 
the controllable SQUID. So, L ■ J <C <fio- On the other 
hand the L ■ J product of the controllable SQUID should 
be sufficiently large, so that the flux expelled is reason- 
ably large. This is particularly important because the 
screening current is in the order of or smaller than the 
smallest critical current in the system. In a 7r-state the 
critical current of the junction is smaller than its maxi- 
mum value in the normal state (see Fig.|J). If these two 
constraints are met we expect a saw tooth behavior of 
the circulating current in the controllable SQUID as a 
function of 4> ex t) a s shown in Fig.|l|. Hence we expect a 
saw tooth behavior of the magnetic field induced in the 
measurement SQUID. 

The device is made on thermally oxidized silicon. In 
the first two steps the measurement SQUID is made us- 
ing conventional e-beam lithography with a double layer 
of PMMA and lift-off. We start with the normal regions 
of the junctions, which are made of 40 nm of gold un- 
der which we use 5 nm of titanium as an adhesion layer. 
In the second step we DC-sputter deposit the niobium 
loop (60 nm) and contacts, which is done after a short 
in situ argon etching of the gold electrodes. The sepa- 
ration between the Nb electrode is 450 nm. To isolate 
the measurement SQUID we cover the entire structure, 
except from the contacts (outside the range of Fig. ||) 
with 170 nm AI2O3 using RF sputtering. In the last two 
steps we fabricate a controllable SQUID on top of the 
small SQUID using the same fabrication procedure as 
described in section III. 

The experiment is done at 1.4 K. The first step is 
to determine the maximum flux generated by the mea- 
surement SQUID. To be able to do this we measure 
the critical current of the the measurement SQUID (at 
4>ext — l/20o and 0) while suppressing the flux gener- 
ated by the controllable SQUID. The latter is achieved 
by means of applying a control voltage exactly equal to 
the critical control voltage to the bottom junction of the 
controllable SQUID. This results in a total suppression of 
the critical current of this junction, i.e. J=0. Hence we 
obtain the critical current of the measurement SQUID at 
(j) ext = of 4.2 \xA and 200 nA at (j) ext = O.50 o . The 
maximum flux generated by the screening current in the 
measurement SQUID can be estimated from these mea- 
surements using Eq. |l|, || and || to be less than O.O50o> 
assuming the worst case scenario of exactly identical crit- 
ical currents for both junctionsES. This indicates that the 
measurement SQUID does not influence the controllable 
SQUID. In the next step we estimate the self inductance 
of the controllable SQUID using the data given in Table 

Fn. 

The critical currents of the bottom and top junctions, 
J C 2 and I c i respectively, are measured by applying pre- 
cisely the critical control voltage to one junction (reduc- 
ing its critical current to 0) and measure the transport 
critical current through the SQUID. From these values 
we can estimate L=65 pH. This value of L would result 
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FIG. 9: Scanning electron microscope of a device that consists 
of a controllable 7T-SQUID with a smaller conventional DC 
SQUID (measurement SQUID) placed inside its loop. The 
two devices are electrically isolated by means of 170 nm of 
AI2O3. The flux expelled from the big controllable SQUID is 
detected using the measurement SQUID. 



TABLE III: Parameters of the controllable SQUID to es- 
timate the self inductance taken at 1.4 K. I c ,max = Ic at 
4>ext — and I c ,min=Ic at <f) ex t = \4>o- Ic2 (lei) is the criti- 
cal current of the bottom (top) junction in Fig|. 



Ic,MAX [fJ.A] 


Ic,MIN [fiA] 


IcB \pA] 


IcT \pA] 


88 


76 


74 
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SQUID over several fluxquanta within less than half a 
fluxquantum in the measurement SQUID. To make the 
result more clear we can extract the magnetic field in- 
duced in the measurement SQUID by the screening cur- 
rent in the controllable SQUID from this data: We sub- 
tract the grey curve from the black one, and multiply 
this difference by (fj^) -1 ; which is obtained by numeri- 
cally differentiating the grey curve. Obviously the data 
around 0, 0o are not obtainable because the 

sensitivity of the measurement SQUID is zero at these 
points. We perform this analysis on a set of data at dif- 
ferent values of V c b ■ The. result is shown in Figjll] for 
2 different values of V^Kll. We observe that the black 
line, which is the same measurement as shown in Fig. 
[Io|(V C 2 > V r C 2,critica;=0.62 mV), is a saw tooth function 
around integer flux, consistent with the behavior of the 
screening current in a SQUID with a relatively large self 
inductance, as shown in Figjljc. The grey line, for which 
V C 2 > V C 2, C ritica/=l-08 mV, shows the exact same behav- 
ior, with however a shift of half a fux quantum. In this 
case a screening current flows, in either the clockwise or 
counterclockwise direction, at integer external flux. 

These measurements indicate that it is possible, 
around an integer flux quantum, to switch between a 
conventional state without a screening current and a 
state with a (bistable) screening current by applying 
V C 2 > V C 2,critical- This is illustrated more clearly by the 
data in Fig.[l2| In this figure we show a measurement 
of the magnetic field induced in the bottom SQUID as a 
function of the V C 2 (in stead of (f> ex t) for two different val- 
ues of the applied external flux around -1 fluxquantum. 
At V C 2 < V C 2,criticai the field induced in the measure- 
ment SQUID is a horizontal line close to 0, indicating 
that hardly any flux is expelled from the controllable 
SQUID. Hence no screening current flows around the 



in a maximum expelled flux at 1.4 K, if no control volt- 
ages are applied to the junctions, of 0.34 </>o associated 
with J = 11/iA 

In the experiment we measure the Vsquid — B oscil- 
lations of the measurement SQUID using a suitable low 
frequency AC bias current. This is done for different val- 
ues of V C 2 on the controllable SQUID. A representative 
result is shown in Fig. [l^. The grey line indicates the 
situation where V C 2 = V C 2 tC ritical = 0.745toV". In this 
case I C 2 — J=Q. The result is an (almost) sinusoidal 
voltage with a period of 0.134 mT, which represents 
cj)Q for the small measurement SQUID. The black line 
shows the situation where V c i < V C 2 critical— 0-62 mV. In 
this case periodic deviations from the grey line are ob- 
served associated with the screening current in control- 
lable SQUID. The periodicity of the deviations is given 
by B((f> . controUab ieSQUiD) = -B($ ) =0-029 mT. The ra- 



tio 



B($o) _ 0.029 



i| which is the ratio of the areas of 



B(4> Q ) 0.134 

the two SQUIDS. It is this difference in area and hence 
in periodicity, that makes it possible to obtain informa- 
tion over the circulating current in the large controllable 
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FIG. 10: Voltage measured over the measurement SQUID as a 
function of magnetic field for two values of the control voltage 
on the bottom junction of the controllable SQUID. V C 2=0.745 
represents the critical control voltage. 
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FIG. 11: The magnetic field induced in the measurement 
SQUID by the flux generated by the screening current in the 
controllable SQUID as a function of the external flux given 
in units of $o, representing one flux quantum through the 
controllable SQUID. The data is given for two different val- 
ues of the control voltage over the bottom junction of the 
controllable SQUID. 



loop of the controllable SQUID as expected. However, 
if V C 2 > V C 2,critical we see an increase (or decrease) of 
the induced magnetic field, associated with the screen- 
ing current caused by the introduction of an extra phase 
factor of 7r (-7r) in the SQUID loop. The sign of the 
signal depends on the exact value of the external flux. 
At 4> ext — —1.100 (black line), the induced field is neg- 
ative, consistent with Fig.[ll]. This clearly indicates that 
it is possible to switch from a state with no current run- 
ning around the SQUID loop to a non zero current state. 
The small a-symmctrics in the figure are caused by the 
magnetic field associated with the current flowing in the 
control channel of the bottom junction. 



VII. CONCLUSIONS 

We have demonstrated that in a controllable 7T-SQUID 
a screening current can be generated without the appli- 
cation of an external magnetic field by putting one con- 



trollable 7r-junction of the SQUID into a 7r-state. At a 
measurement temperature of 1.4 K and with the present 
geometry, only a Nb-Ag controllable tt junction shows 
the transition to a 1 state, in contrast with a Nb-Cu 
device. We have shown that this difference is caused 
by the difference in electron-electron interaction strength 
between Ag and Cu, consistent with direct measurements 
of this quantity on similar diffusive wires using tunnel 
junction spectroscopyEJEII. The possibility of adding an 
extra phase of it to the SQUID loop is demonstrated not 
only by the screening current, but also by the fact that 
the Vsquid—B oscillations can be shifted by a phase fac- 
tor of 7r or 2tt, depending wether one or two controllable 
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FIG. 12: The magnetic field induced in the measurement 
SQUID by the flux generated by the screening current in 
the controllable SQUID as a function of the control voltage 
applied over the bottom control channel of the controllable 
SQUID. The data is shown for 2 values of the external flux $ 
through the controllable SQUID around -1 fluxquantum. 



junctions are put on a 7r-state. 

We gratefully acknowledge H. Pothier and F.K. Wil- 
helm for discussions and for making their computer pro- 
grams available to us. This work was supported by the 
Nederlandse Organisatie voor Wetenschappelijk Onder- 
zoek (NWO) through the Stichting voor Fundamenteel 
Onderzoek der Materie (FOM). 



1 A.F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). 

2 A.F. Volkov, H. Takayanagi, Phys. Rev. B. 56, 11184 
(1997). 

3 F.K. Wilhelm, G. Schon, A.D. Zaikin, Phys. Rev. Lett. 81, 
1682 (1998). 

4 S.-K Yip, Phys. Rev. B 58, 5803 (1998). 

5 B.J. van Wees, K.-M.H.Lenssen and C.J.P.M.Harmans 
Phys. Rev. B Rap. Com. 44, 470 1991 

6 A.F. Morpurgo, T.M. Klapwijk and B.J. van Wees, Appl. 
Phys. Lett. 72, 966-968 (1998) 



7 K. Neurohr, Th. Schapers, J. Malindretos, S. Lachenmann, 
A.I. Braginski, H. Luth, M. Behet, G. Borghs, A. A. Gol- 
ubov, Phys. Rev. B 59, 11197 (1999). 

8 J. Kutchinsky, R. Taboryski, C.B. Sorensen, J.B. Hansen 
and P.E. Lindelof, Phys. Rev. Let. 83, 4856 (1999) 

9 J. J. A. Baselmans, A.F. Morpurgo, B.J van Wees, T.M. 
Klapwijk, Nature 397, 43 (1999). 

R.Shaikhaidarov, A.F. Volkov, H. Takayanagi, 

V.TPetrashov, P.Delsing, Phys.Rev. B 62. R14649 
(2000) 



11 



C. C. Tsuei, J.R. Kirtley, C.C. Chi, Lock See Yu-Jahnes, 
A. Gupta, T. Shaw, J.Z. SUn, M.B. Ketchen, Phys. Rev. 
Let. 73, 593 (1994); A. Mathai, Y. Gim, R.C. Black, A. 
Amar, F.C. Wellstood, Phys. Rev. Let. 74, 4523 (1995) 

D. J. van Harlingen, Rev. Mod. Phys. 67,515 (1995) 
V.V. Ryazanov, V.A. Oboznov, A. Yu. Rusanov, A.V. 
Veretennikov, A. A. Golubov, J.Aarts, Phys. Rev. Lett 
86,2427 (2001) 

A system with a weak link between two Superfluid 3He 
reservoirs shows a stable state around a phase <f>=Q and 
a meta stable state around <j) = n with respect to the 
mass current through the weak link. See S.Backhaus, S. 
Pereverzev, R.W. Simmonds, A. Loshak, J.C. Davids and 
R.E Packard, Nature 392, 687 (1998) 
J.Clarke in SQUID sensors, Fundamentals, Fabrication, 
and Applications, editted by H. Weinstock, NATO ASI Se- 
ries 329, Dordrecht, 1-62 (1996). 

M. Cryot and D. Pavuna, in Introduction to superconduc- 
tivity and high-T—c materials , World Scientific, Singa- 
pore, (1992). 

M. Tinkam, in Introduction to Superconductivity, McGraw- 
Hill, Singapore, (1996). 

We have used the fact that the critical current generally 
equals the supercurrent if the phase difference is n/2. This 
is not fully correotc-ibut only leads to minor changes in the 
calculated curvesBu. 

P. Dubos, H. Courtois, B. Pannetier, F.K Wilhelm, A.D. 
Zaikin, and G. Schon, Phys. Rev. B, 63, 064502 (2001) 
H. Pothier, S.Gueron, N.O. Birge, D. Esteve, M.H. De- 
voret, Phys. Rev. Let. 79, 3490 (1997). 
F. Pierre, H. Pothier, D. Esteve, M.H. Devoret, J. of low 
Temp. Phys. 118, 437 (2000). 

A. B. Gougam, F. Pierre, H. Pothier, D. Esteve, Norman 
O. Birge, M.H. Devoret, J. of Low Temp. Phys. 118, 447 
(2000). 

F. Pierre, H. Pothier, D. Esteve, M.H. Devoret, A.B. 
Gougam and Norman O. Birge in Kondo effect and de- 
phasing in low dimensional metallic systems, eds. Venkat 
Chandrasekhar, Chris Van Haesendonck, and Alfred Za- 
wadowski fKlu wer Academic Publi shers. Dordrecht, 2001). 
e-print: arXiv: bond-mat/0012038| (2000). 

B. L. Altshuler, A.G. Aronov, Zh. Eksp. Teor. Fiz 75, 1610 
(1978) [Sov. Phys. JETP 48,812 (1978)] 

B.L. Altshuler A.G. Aronov in Electron- Electron interac- 
tions in Disordered Systems, eds. A.L. Efros, M. Pollak, 
Elsevier Science publishers B.V. (North-Holland, Amster- 
dam 1985) 

F. Pierre in Interactions electron-electron dans les fils 
mesoscopques (electron-electron interactions in mesoscopic 
wires) Ph.D. Thesis CEA-SACLAY (2000) 



F. C. Wellstood, C. Urbina, J. Clarke, Phys. Rev. B 49, 
5942 (1994). 

M. Henny, S. Oberholzer, C. Strunk , C. Schonenberger, 
Phys. Rev. B 59, 2871 (1999) 

The dominant thermal resistance in our system is the one 
associated with electron heat diffusion from the control 
channel into the reservoirsE3. We estimated the electron 
diffusion heating parameter b valid for our samples at 1.4 
K, b < 0.05fc/mV. The calculated effective electron tem- 
perature in the reservoirs T e ff ■ 



yields a 
used in the 



negligible electron heating in the range of V c 
experiments < V c < 1.2 mV. 

J. J. A. Baselmans, B.J. van Wees, T.M. Klapwijk, Phys. 
Rev. B 63, 094504 (2000) 

The use of a Ge suspended mask, which is more stable than 
the presented PMMA mask, is not suitable because the 
electron beam exposure of the large area reservoirs has to 
be done using relatively high beam currents ( nA), to come 
to a reasonable exposure time. This results in stress in the 
PMMA-MA-Ge-PMMA resist, causing the germanium to 
crack. 

D. Esteve et al. in Mesoscopic electron transport: NATO 

ASI series E, vol. 345, 469 (eds. Sohn, L.L. et al.)(Kluwer 

Academic, Dordrecht/Boston/London, 1997) 

B.J. van Wees, H. Takayanagi, in Mesoscopic electron 

transport: NATO ASI series E, vol. 345, 375 (eds. Sohn, 

L.L. et al.)(Kluwer Academic, Dordrecht/Boston/London, 

1997) 

We calculated I c as a function of td with identical param- 
eters and found that at rn=2ns the trasnition to a 7r-state 
would occur at 1.4K. At T < IK the transition would even 
occur with the sample as it is. 

J. J. A. Baselmans, B.J. van Wees, T.M. Klapwijk, Appl. 
Phys. Lett 79, 2940 (2001) 

The self inductance would be 45 pH in the case of identical 
junctions. However, the difference in transport critical cur- 
rent between (j> — and cj> — 4>o/2 can also be attributed to 
small difference in critical currents between the junctions. 
This assumption gives L=0 and a difference in critical cur- 
rent of 10 percent. We must conclude that L<45pH for the 
measurement SQUID. 

In Fig. [ll] we have corrected for a slight horizontal and ver- 
tical shift between the curves that is caused by the current 
associated with V C 2 flowing through the control channel. 
The flux associated with this current causes a horizontal 
shift between the curves due to the coupling in the control- 
lable SQUID, identical to the one in Fig.fr] and a vertical 
shift due to the coupling into the measurement SQUID. 



